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Q ■ Abstract 

Inspired by Quantum Mechanics, we reformulate Hilbert's tenth problem in the 
domain of integer arithmetics into either a problem involving a set of infinitely cou- 
pled differential equations or a problem involving a Shrodinger propagator with some 
, appropriate kernel. Either way, Mathematics and Physics could be combined for 

I I Hilbert's tenth problem and for the notion of effective computability. 

1 Introduction 

The twentieth century witnessed the remarkable discovery of the limits of Mathemat- 
^ I ics, established within itself, through the noncomputable/undecidable results of Hilbert's 
■ tenth problem, Godel's incompleteness theorem, Turing halting problem, and their various 
extensions (see, for example, [1, 2]). Such noncomput ability and undecidability set the 
boundary for computation carried out by mechanical processes, and in doing so it help us 
to understand much better what can be so computed mathematically. 

We have proposed elsewhere [3] a quantum algorithmic approach for the non-computable 
Hilbert's tenth problem [1, 4], which is equivalent to the Turing halting problem and in- 
timately links to the concept of effective computability as defined by the Church- Turing 
thesis. While the proposal is about some quantum processes to be implemented physically, 
it illustrates the surprisingly important role of Physics in the study of computability. This 
is an unusual state of affairs when Physics, which has its roots in the physical world out 
there, could perhaps help setting the limits of Mathematics. 

Inspired by this connection between the two, we present in this work some mathematical 
reformulation of Hilbert's tenth problem. The reformulation is made possible since physical 
theories in general, and Quantum Mechanics in particular, have enjoyed the support and 
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rigour of mathematical languages. We wish to stress here that, despite of the inspiration, 
the connection is established entirely in the domain of Mathematics; we need not appeal 
to some real physical processes as we do with the proposed quantum algorithms in [3]. It 
is hoped that such reformulation may lead to new insights of the problem. 

In the next section we briefly state the problem and the inspired connection with 
operators acting on some infinite-dimensional Hilbert space. ^Prom this we first derive a set 
of non-lincarly, coupled differential equations (eqs. (13, 14)) then also a linear Scrhodinger 
equation (eq. (30)), each of which cases separately contains the sought-after decision result. 
If one could find a universal procedure to derive certain information from these differential 
equations (not necessarily by solving them explicitly but could be by other means) then 
one would have settled Hilbert's tenth problem in the positive. 

2 Hilbert's tenth problem and Hilbert space 

At the turn of the last century, David Hilbert listed 25 important problems, among which 
the problem number ten could be rephrased as: 

Given any polynomial equation with any number of unknowns and with in- 
teger coefficients: To devise a universal process according to which it can be 
determined by a finite number of operations whether the equation has integer 
solutions. 

This decision problem for such polynomial equations, which are also known as Diophan- 
tine equations, has eventually been shown in 1970 by Matiyascvich to be undecidable [4, 1] 
in the Turing sense. It is consequently noncomputable/undecidable in the most general 
sense if one accepts, as almost everyone does, the Church- Turing thesis of comput ability. 
Since exponential Diophantine, with the unknowns in the exponents as in the example 
of the Fermat's last theorem, can be shown to be Diophantine with supplementary equa- 
tions, the study of Diophantine equations essentially covers the class of partial recursive 
functions, which is at the foundation of computability. The undecidability result is thus 
singularly important: Hilbert's tenth problem could be solved if and only if could be the 
Turing halting problem. (See [5] also.) 

Given a Diophantine equation with K unknowns x's 



it suffices in general to consider the existence of nonnegative integer solutions. 

Following [3] we link the equation above with the following hermitean operator acting 
on some appropriate Fock space (a special type of Hilbert space) 




(1) 




(2) 



where 



K, %] = 0. 
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The Fock space is built out of the "vacuum" |0a) by repeating apphcations of the creation 

operators a]. 

The operator (2) has a semidcfinitc and discrete spectrum (D(ni, ■ ■ • , n^))^. This 
spectrum has an eigenstate \Eg) corresponding to the smaUest eigenvalue Eg. (If the 
hermitean operator is considered as a hamiltonian for some dynamical process then these 
are the ground state and its energy, respectively.) 

It is clear that the Diophantine equation (1) has at least one integer solution if and 
only i{Eg = Q. 

To sort out this Eg among the infinitely many eigenvalues is almost an impossible task. 
The trick we will use, as inspired by quantum adiabatic processes, is to tag the state \Eg) 
by some other known state \Ei) which is the ground state of some other operator Hj and 
can be smoothly connected to \Eg) through some continuous parameter s e [0, 1]. To that 
end, we consider the interpolating operator 

rL{s) = Hi + f{s)iHp-Hi), 

= Hi + f{s)W, (3) 

which has an eigenproblem at each instant s, 

[n{s)-Eg{s)]\E,{s))^0, g = 0,l,--. (4) 

with the subscript ordering of the sizes of the eigenvalues, and /(s) some continuous and 
monotonically increasing function in [0, 1] 

/(O) = 0; /(I) = 1. (5) 

Clearly, £"0(0) = Ej and Eo{l) = Eg. It turns out that for the function Eo{s) to connect a 
ground state to another ground state we require that 

[Hp, Hi] 0, (6) 

and that Hj has a discrete spectrum with a non-degenerate ground state \Ei). An example 
of such Hi is, which we will employ from now on, 

K 

Hi ^^{al-a*){ai-ai), (7) 

i=l 

in which case, Ei = and \Ei) = \ai - ■ ■ ax), that is, the coherent state 

l«) = e-— (8) 
n=o ynl 

where a e C and \na) are the eigenstates of a^a with eigenvalues n. 
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3 The spectral flow 

We now derive the differential equations for tlie tagging connection for tlie instantaneous 
eigenvalues and eigenvectors at different instant s in (4). 

Note firstly that, from the normalisation condition {Eq\Eq) — 1, we can write 

{Eg\ds\Eg) = -ids(j)q, 

for some real (pg. This can be absorbed away with the redefinition 

e^^^W|E,(.)) - \Eg{s)), 

upon which 

{Eg\ds\Eg) = 0. (9) 

(The phases (pg are related to Berry's phase in a dynamical process.) 
Differentiating (4) with respect to s yields 

[f'{s)W-d,Eg]\Eg)+[n-Eg]d,\Eg) = 0. (lO) 

We next insert the resolution of unity at each instant s, 

oo 

1 = ^\E^{s)){E^{s)\, 

m=0 

just after 7i in (10) to get, vy virture of (9), 

oo 

Egd,\Eg) = [f'{s)W-d,Eg]\Eg)+j2E^{Ejds\Eg)\Ej. (ll) 

The inner product of the last equation with \Ei) gives 

{Eg-El){El\d,\Eg) = f'{s){El\W\Eg)-dsEgdgl. (l2) 

Thus, for g 7^ / this gives the components of ds\Eq) in \Ei), provided Eg ^ Ei at any 
s G (0,1), a condition we will investigate in the next section. Consequently, together 
with (9), 

a.\E,) = ns)t^-f^\E,). (13) 

l^g Eg -El 



Also, putting g = Hn (12) we have 

dsEg{s) = f{s){Eg{s)\W\Eg{s)). (l4) 

Equations (13) and (14) form the set of infinitely coupled differential equations providing 
the tagging linkage we have been after. 
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Analytical and numerical methods could now be employed to investigate the unknown 
ground state of Hp from the constructively known spectrum of Hj as the initial conditions. 
In this reformulation, the Diophantine equation (1) has at least one integer solution if and 
only if 

lim£;o(s) = 0. (15) 

The limiting process is necessary since Hp, i.e. ?i(l), will have a degenerate spectrum 
because of certain symmetry [Hp commutes with a|aj). 

The equations above are infinitely coupled and cannot be solved explicitly in general. 
But we are only interested in certain information about the ground state. And since the 
influence on the ground state by states having larger and larger indices diminishes more and 
more thanks to the denominators in (13) (once no degeneracy is assured), this information 
may be derived, numerically or otherwise, with some truncation to a finite number of states 
involved. The size of the truncation cannot be universal and is of course dependent on the 
particular Diophantine equation under consideration. 



4 Expansion in the number basis 



We derive in this section an explicit form for the equations (13, 14) in the case of 2 
unknowns, i.e. i^' = 2 in (1). 

At any instant s we expand the eigenvectors \Eq[s)) in the basis consisting of the states 
\n1n2) of the number operators a\ai and 020.2, 



\Eq{s)) = ^ Cg;nin2(s)|nin2), 
rai,"2=0 



(16) 



upon which all the s-dependency are now carried by the coefficient functions Cq-^nin2 (s) of 
the expansion. Direct substitution of the expansion into (13, 14) gives 

dsEg{s) = f'{s) J2 {|C<?;nin2(s)P D^{ni,n2) - rii - n2 - |ai|^ - jaal 



and 



9s^q;nin2 (^) 



+ aWni + lC*.(„^+i)„^(s)Cq;„,„2(s) + c.c. 
+ a2Vn2 + lC;ni(n2+i)(«)C9;nin2(s) + cc} ; 

{Q*nin2('5)C<7;nin2(s) D'^{ni,n2) - Tli - ^2 - - |a2| 

(ni+l)n2 {^) 

+ a2^Jn2 + lC,*„^(„2+i)(s)Cg;„i„2(s) 

+ ayn2 + lC/*„i„2(s)Cq;ni(n2+l)(s)} 



(17) 



(18) 
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The appropriate initial conditions for these infinitely coupled differential equations can 
be derived from the eigenvalues and eigenvectors of the initial hermitean operator Hj. If 
we choose to index the initial eigenvalues as 

{Eo{0),E,{0),E2{0),Es{0),---) = (0,1,1,2,-.-), (19) 

then the coefficient functions at s = can be inferred as follows. ^Prom the expression for 
the coherent state (8) 



\Eo{0)) = \a1a2), 



„ "1 „ "2 
ai a2 



= Yl ^==l^l^2), (20) 

ni,n2=0 VnM 



we arrive at 

Co;n.nM = e-^^^L=4=. (21) 

The next excited states arc doubly degenerate, Ei{Q) = £"2(0) = 1, 

\Ei{0)) = {a\-al)\aia2), 

= -allEom + e-''^ E ^^V^l\{n, + l)n2); (22) 

ni,n2=o Vnilna! 



and 

[£'2(0)) = {al - a2)\aia2) , 



= -a;|£o(0)) + e-^^ E ^^^^A^|nl(r^2 + l)). (23) 



ni,n2 

The last two equations subsequently yield 



-Q;tC'o;on2(0); ni = 0. 



And 



<^2;nin2(0) = < V"i'("2-1)! (25) 

I -Qi2Co;mo(0); 712 = 0. 

The expressions (19, 21, 24, 25), etc . . . are the initial conditions at s = for the differ- 
ential equations (17, 18). However, because of the degeneracy in the initial conditions (19) 
we may have to integrate the differential equations from some s = e infinitesimally away 
from zero - in which case the degeneracy should be lifted (as justified in the next section, 
there is no generic level crossing in the open interval < s < 1); and where the new initial 
conditions can be estimated numerically (to be very closed to the values at s = 0) by the 
familiar perturbation theory in quantum mechanics. 
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5 No crossing for a single pair of levels 



We now argue that at any instant sq e (0, 1) if there is only one pair of eigenvalues, say 
(1,1+1) out of the infinitely many, can come very close together then they can never actually 
cross. The arguments are similar to those of perturbation theory for nearly degenerate 
levels [6]. 

As can be seen from (13), in a neighbourhood around Sq we need only consider the two 
states \Ei) and \Ei+i) as they are so strongly coupled that the rest can be safely ignored. 
As our world just becomes a two-dimensional space, we can linearly decompose the two 
states at the next instant (so + f^s) in terms of the two at Sq. Now with this decomposition, 
we need only to solve a two-by-two determinant for the eigenvalue problem for Ti.{so + 5s). 
The end result for the gap is 

^i,i+i{so + Ss) = Ei+i{so + Ss) - Ei{so + Ss), 

= ^[A;,;+i(5o) + dsf'{so){Wi+^,i+^ - Wi^i)]-^ + 4Ss^\f'{so)Wi,i+^\^, 

(26) 

where Wij = {Ei{so)\W\Ej{so)) . The matrix element in the last term in the square root 
is just proportional to the difference 

\{Ei{so)\niso + 5s)-niso)\Ei+,{so))\^ = \{Ei{so)\n{so + 5s)\Ei+iiso))\^ . (27) 

Since 

[n{so + Ss),n{so)] = 5sf{so)[Hp,Hj]+0{Ss'')^0, (28) 

even restricting to our two-dimensional subspace. Then, as the eigenvectors of Tt{so) cannot 
diagonalise Tl{so + 5s), upon which the off diagonal elements (27) cannot be zero. Thus the 
gap (26) can never be zero, no matter how vanishingly small it is at the previous instant. 
(If the commutator (28) vanishes for our two-dimensional matrices, in contradiction to its 
nonvanishing for the operators in the whole Hilbert space, then our assumption of being 
able to ignore all other eigenstates apart from the two in consideration is not valid-perhaps 
because the gap at sq is not yet small enough.) 

We note that from (26) or (14) we can derive a differential equation for the gap and 
look for its minimum value. 

The above arguments are only applicable for "accidental degeneracy", that is when 
there is exactly one single pair of levels comes very closed. The reasoning fails when we are 
not allowed to isolate a two-dimensional subspace and treat it separately as in the above. 
Namely, it may fail when there are not two but three or more levels come crossing at one 
point; then the denominators in (13) force us to consider a larger subspace. The resulted 
determinant will have larger dimensions and not enough constraints to keep all the gaps 
non-zero. 

The above arguments for level avoidance also fail when there is another crossing pair, 
say {q, q + 1), elsewhere in the spectrum at that same instant. Then the feedback of that 
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pair {q, q + 1) through their states to the right hand side of (13) for the pair {I, I + 1) may 
not be ignored; and we end up with a dimensionally larger subspace again. 

But it is neither accidental nor typical that three or more levels cross at exactly one 
point, or two or more pairs become degenerate at exactly the same instant. Those events 
belong to a zero-measure set of events, unless there must be a reason. It should be a sym- 
metry reason, that is, 7i{so) should commute with some other hermitean operator (s). Our 
mathematical elaboration above agrees with the observation of symmetry and degeneracy 
in [8]. 

With care we can slightly modify the derivation for (13, 14) to come up with similar 
equations even when there is some degeneracy in [0, 1]. But for the condition (15) to be 
the indicator for the existence of solutions of the Diophantine equation, simple topological 
consideration only requires that the initial ground state |£'o(0)) = \Ej) is not degenerate 
and that this state doesnot cross with any other state in the open interval s G (0, 1). With 
the freedom of choice for Hj satisfying (6), we should be able to eliminate any symmetry 
in the open interval s e (0, 1) for H{s) in order to have a stronger condition of totally 
avoided crossing. This is because that accidental symmetry cannot persist with widely 
different choices of a,j in different starting HjS (which result in different sets of differential 
equations that can be used to cross confirm each other). 

Alternatively, one could also systematically eliminate the degeneracy effects, if any, 
caused by some accidental symmetry for a general Diophantine equation by considering a 
modified Hp 

Hp ^ //p + ^(e.aj + e*a,), (29) 

i=l 

with |ej| <^ 1. In the limits \ei\ — > 0, we recover our original Hp and would also be able to, 
by going through the limiting processes, discover and eliminate any accidental symmetry. 
It is clearly seen that the extra terms in (29) will also remove any degeneracy of 7Y(1) = Hp. 
(This is similar to well-exploited physical trick in atomic physics of degeneracy removing 
by small biased external fields.) 

6 Hilbert's tenth and the Schrodinger equation 

The decision result for Hilbert's tenth problem is also contained in yet another type of 
differential equation, apart from the nonlinear equations (13, 14). The linear equation is 
just the Schrodinger equation which captures the dynamics of our quantum algorithm [3]. 
Let \ip{t)) be the quantum state at time t, its time evolution is given in quantum mechanics 
by the equation 

dt\m) = -tn{t/TMit)), (30) 

1^(0)) = \ai---aK), 
where we have chosen the initial state at time i = to be the ground state of Hj. 
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Unlike the case of nonlinear differential equations of a previous section, now we could 
try to make use of a powerful computability result in analysis which is known as the First 
Main Theorem [9]. Essentially, the Theorem states that a bounded linear operator from a 
Banach space to a Banach space which maps a computable sequence of spanning vectors 
into another computable sequence will also map any computable element into another 
computable element. For the case at hand, our Schrodinger equation defines a linear 
operator. 



where T is the time-ordering symbol, mapping the initial state to the final state in the same 
separable Hilbert space. Now, our initial state \ai • • - ax) is computable by construction. 
On the other hand, the linear operator (31) coming out of the Schrodinger equation should 
be unitary and thus be bounded. Hence, the only remaining condition of the Theorem 
to be checked is whether the image of a particular computable basis is computable or 
not. Speculatively, if it is the case then Hilbert 's tenth problem of integer arithmetic 
is dccidablc through the use of mathematical analysis tools (plus those of the theory of 
infinite-dimensional operators)! This will be investigated elsewhere. 

Nonetheless, we speculate that for the Schrodinger equation to offer some new results 
here then there must be no level crossing in the spectral flow. But here as we start with the 
ground state of Hi we only require no level crossing for the instantaneous ground state, 
unlike the situation with the nonlinear equations previously where we have required no 
crossing for all levels. Adapting Ruskai's arguments [10] we can show that the ground 
state of (3) is non-degenerate for s e (0, 1). (The arguments are only applicable for the 
ground state but, interestingly, the conclusion of the last section can also be supported by 
them. The trick is to use Ruskai's arguments for the two-dimensional "subspace" there; as 
there are only two levels, the lower level is now the ground state. Once again, this trick 
is not applicable when we cannot, because of some symmetry reason as already discussed 
in the last section, isolate such two-dimensional subspaces. With more dimensions than 
two, adaptation of Ruskai's cannot rule out the crossing of other states different from the 
ground state.) 

It should also be noted that the above quantum mechanical approach to Turing- 
noncomputable problems is in contrast to the claim in [11] that quantum Turing machines 
compute exactly the same class of functions, albeit perhaps more efficiently, which can be 
computed by classical Turing machines. However, the quantum Turing machine approach 
is a direct generalisation of that of the classical Turing machines but with qubits and some 
universal set of one-qubit and two-qubit unitary gates to build up, step by step, dimension- 
ally larger, but still dimensionally finite unitary operations. This universal set is chosen 
on its ability to evaluate any desirable classical logic function. Our approach, on the other 
hand, is from the start based on infinite-dimension Hamiltonians acting on some Fock 
space and also based on the special properties and unique status of their ground states. 
The unitary operations are then the Schrodinger time evolutions. The infinite dimension- 
ality together with the unique energetic status of the vacuum could be the reasons behind 



U{T) 




(31) 
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the ability to compute, in a finite number of steps, what the dimensionally finite unitary 
operators of the standard quantum Turing computation cannot do in a finite number of 
steps. Note that it was the general Hamiltonian computation that was discussed by Benioff 
and Feynman [12, 13] in the conception days of quantum computation. 

Indeed, Nielsen [14] has also found no logical contradiction in applying the most general 
quantum mechanical principles to the computation of the classical noncomputable, unless 
certain Hermitean operators cannot somehow be realised as observables or certain unitary 
processes cannot somehow be admitted as quantum dynamics. And up to now we do not 
have any evidence nor any principles that prohibit these kinds of observables and dynamics. 

7 Concluding remarks 

Inspired by Quantum Mechanics, we have reformulated the question of solution existence 
of a Diophantine equation into the question of certain properties contained in an infinitely 
coupled set of differential equations. In words, we encode the answer of the former question 
into the smallest eigenvalue and corresponding eigenvector of a hermitean operator whose 
integer- valued spectrum is bounded from below. And to find these eigen-properties we next 
deform the operator continuously to another operator whose spectrum is known. Once 
the deformation is also expressible in the form of a set of nonlinearly coupled differential 
equations, we could now start from the constructive knowns as a handle to study the 
desired unknowns. 

In addition, we also explicitly present a linear Schrodinger equation whose solution 
at some time T from an appropriate initial condition contains the information about the 
decision result for the Diophantine equation under investigation. 

Note that these reformulations are entirely based on mathematics. If a general mathe- 
matical method could be found to extract the required information from these differential 
equations for any given Diophantine equation then one would have the solution to Hilbert's 
tenth problem itself. This may be unlikely but not be as contradictory as it seems -because 
the unsolvability of Hilbert's tenth problem is only established in the framework of inte- 
ger arithmetic and in Turing computability, not necessarily in Mathematics in general. 
Tarski [7] has shown that the question about the existence of real solutions of polynomials 
over the reals is, in fact, decidable. 

In the case of the linear Schrodinger (30), we could also exploit a powerful computabil- 
ity result in analysis which is known as the First Main Theorem [9]. With this Theorem, 
computability has been illustrated to be indeed context/framework-dependent in the ex- 
ample of the classical wave equation: whether the solution to this equation is computable 
or not depends crucially on the initial functions and definitions of the norm employed. 
We refer to the original literature for a thorough discussion of this remarkable property. 
In order to establish the mathematical decidability of Hilbert's problem in this particular 
reformulation, we will need to investigate (and this will be done elsewhere) whether the 
unitary transformation (31) satisfies all the conditions of the Theorem or not. 

On the other hand, even though this Theorem for linear operators is not applicable to 
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the set of non-linear equations (13, 14), such reformulation of Hilbert's tenth problem with 
continuous variable has opened up many new directions for further investigations. 
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